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We present a brief review of CPV and mixing measurements in the charm sector, 
with emphasys in results published or presented since the previous edition of the 
Physics in Collision Symposia. 

1 Introduction 

Neutral mesons mixing occur in nature due to the fact that the mass eigenstates 
are not weak interaction eigenstates. An initially pure beam of D° or D° mesons 
evolves in time to a mixture of D° and D° states with rates determined by the 
solution of the Schroedinger equation 

i |(^o)(<)=H« ff (go)(*)> (1) 

where the diagonal elements of the effective Hamiltonian H e fr = M — ±T and of 
the hermitian matrices M and T are related to the flavour conserving D°(D°) —> 
D°(D°) transitions, while the non-diagonal ones are related to the flavour changing 
D° <H- D° transitions. The solution of Eqjljare the eigenstates |-Di,2) = p\D°) ± 
q\D°), with complex coefficients peg satisfying \/|f>| 2 + \q\ 2 = 1- The eigenstates 
have well defined masses and widths m,i,Ti and 7712,^2- The parameters which 
govern the mixing are x = ml ~ m2 and y = ^jf 2 -, where T = Fl ^" r2 is the mean 
decay width. For charm mesons, both x and y are small. There are essentially two 
approaches used to predict their values in the Standard Model (SM), none of them 
leading to fully reliable results (see, for example, 

The transitions leading to the D° — D° oscillations are those represented in Fig. [1] 
Short distance contributions are illustrated by the box diagram on the left and 
long distance transitions, involving intermediate states accessible by both D° and 
D°, on the right. As seen from these diagrams, charm mixing gives us access to 
transitions involving intermediate states with down-type quarks. In fact, this is the 
only system which offers this probe, since the top quark is too heavy to hadronize 
and the 7r° is its own anti-particle. 

While x receives contributions from both the short distance and long distance 
processes, y only receives contributions from the long distance transitions. Since 
x is the only one to receive contributions from virtual intermediate states, new 
physics introducing off-shell particles in the box diagrams would only affect x and 
therefore it is a sort of consensus that an unambiguous sign of NP would be given 
by a x value well above y. The first evidences of mixing were obtained in 2007 
by BaBar [4j, Belle [5] and CDF jfjj. At present, both x and y are constrained by 
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1|2|3 and references therein). 
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several measurements at the values x = (0.65+ lg )% and y — (0.73±0.12)% and the 
no mixing hypothesis is excluded with a statistical significance above 10 standard 
deviations (a) [7], under the assumption of no CP violation. Although theoretical 
predictions are not precise, the experimental values of the mixing parameters are 
now consistent with SM expectations [8], even if still considered at the upper level 
of possible values (9j. 
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Figure 1. Short and long distance contributions to D° — D° mixing. 



The CP symmetry can be violated in the mixing, if \q\ ^ \p\; in the decay, if 
the magnitudes of the instantaneous decay amplitudes Af = (f\H\D°) and Af — 
(/|H|-D°) for CP conjugate processes are not the same; and, for decays reachable 
by both the D° and D°, in the interference between mixing and decay. This kind 
of mixing induced CPV can occur even in the absence of CPV in the mixing and 
in the decay, if there is a non null phase difference (f> between the mixing and the 
decay amplitudes. 

In the SM, CP violation (CPV) is suppressed by the small values of the CKM 



parameters V u b and V c b- In the mixing, it is expected to be at the level of 10 -4 [10 
independently of the final state. Direct CPV can only occur in singly Cabibbo 
suppressed (SCS) decays, since their amplitudes can receive contributions from 
processes with virtual b-quarks, involving the complex elements of the CKM matrix. 
Existing predictions for SCS modes in the SM were, at least until recently, of the 
order of 10 -3 or less. The interest in CPV in charm decays is, therefore, related to 
searches for new physics effects in a SM suppressed environment. Experimentally, 
the only existing evidence of CPV in the charm sector is the measurement of the 
difference between the time integrated asymmetries of the decays D° —> K + K~ 
that will be further discussed in this document. 
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A very good review on CPV and mixing on the charm sector is given in [3j. 
Here we report on some recent measurements performed by BaBar, Belle, CDF 
and LHCb since the previous edition of the Physics in Collision conference series. 



2 Time dependent measurements of D — D° mixing and CP 
violation 

There are several ways to access the mixing parameters x and y experimentally, 
using decays to various final states. Each method provide different sensitivities 
to these parameters and, by combining these measurements, we can improve our 
knowledge about them. The HFAG group [12] provides world averages and best fit 
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values to these parameters using results provided by several experiments at different 
machines around the world. A good overview of the most up to date averages are 
given in |7j. 

One of the most sensitive ways of determining x and y is via a time dependent 
amplitude analysis of the decays D° — > K s tt + tt~ and D° — > K S K + K~ . The 
two final states are CP eigenstates reachable by both the D° and D°. The time 
dependent rate for these decays is proportional to terms arising from the interference 
between the direct decay and the decay via mixing: 

K(t) cx |Ai| 2 e^* + \A 2 \ 2 e- yt + 2Re[A 1 A* 2 ] cos(xt) + 2Im\A 1 Al] sin{xt) (2) 

where A\j, — A % A and A and A are the instantaneous decay amplitudes for D° 
and D°, written as a function of the 2-body invariant masses M + (KgK + ) and M~ 
(KgK~). In order to determine the flavour of the meson at production, it is required 
to form a vertex with a low momentum pion tt s , consistent with the decay vertex of 
a D* + -t 7T+D or D*~ -t n~D°. The charge of the ir s identifies the D° flavour. 



This analysis was first performed by CLEO [13] and subsequently by Belle 14 and 



BaBar 15 , with much larger samples. An update of the D° — > K s tt + it~ analysis 
by Belle was presented this summer, using their full data set of 920 fb _1 and an 
improved tracking, which provides a signal yield of 1.23 M events with a purity of 



95.6% 16 . They parametrize the amplitudes A and A as a sum of 12 quasi-two- 
body amplitudes describing the P and D waves and two components to describe 
the Ktt (LASS model) and the nir (K-matrix model) s-waves. From a maximum 
likelihood fit to the data in the signal region, the mixing parameters x and y, the 
average -meson lifetime tjjo and the magnitudes and phases of the resonances are 
directly extracted. The preliminary results shown are x = (0.56 ± 0.19jQ g 03 jQ 9 06 )% 
and y = (0.30 ± 0.15^Q Qgj;g'Qg)%, where the first uncertainty is statistical, the 
second is due to experimental systematic effects and the third systematical due to 
the model used to fit the Dalitz plot. At the time of the conference, this was the 
most precise single measurement of the mixing parameters. 

Another measurement which is sensitive to mixing is the deviation from unity of 
the ratio of the effective lifetimes measured in decays to the CP eigenstates K~K + 
and 7r - 7r + with respect to the CP mixed state K~ir + . Since mixing is very slow, 
the decay time distributions are, to a good approximation, simply exponential. The 
fits to the decay time distributions of the different samples use a model where signal 
events are described by a single exponential and provide the effective lifetimes used 
to build the observable 

VCP = t(D° ^K- k J) ~ 1 ~ 5 + I « \) ycos ^~ (IpI " If I) xsin< ^ W 

where h~h + corresponds to 7r _ 7r + or K~ K + , the phase <f> is considered universal 
and a contribution from CPV in the decay is neglected. In the absence of CPV, 
\q/p\ — 1 and (f> — 0, hence yep = y- This measurement was previously performed 
by several collaborations and provided the first evidence for mixing obtained by the 
BELLE experiment in 2007 IB]. In April 2011, both BaBar and Belle presented new 
results using their full data samples, with integrated luminosities of 468 fb -1 17 
and 976 fb -1 [l8], respectively. LHCb also published a first result for yep us 



ing 29 pb of data collected in 2010 19 . In addition to the D* ± tagged samples, 
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BaBar used also untagged samples of D°(D°) -> if ± vr =F and D°(D°) -> 
decays. These results are presented in Table[T] together with the size of the samples 
used in each measurement. 

The effective lifetimes obtained by fitting the decay time distributions of the 
tagged samples separately for the CP conjugate decays provide a second quantity, 
Ar, sensitive to CPV: 

Ar(h~k + ) ~ T{D ° -> h ~ h+) - T{D ° k ' h+) ~ 1 \(\i\-\ZWv COS <b-m + \*\\xsm 
A T {h h )- T(D o^ h - h+)+T{D o^ h - h+) ~2 LCI ^ I \<,\)y™s<P [\,\ + \ q \)**» 

(4) 

considering as valid the same approximations used for yep- 

The three results for yep and Ar are consistent with each other and with 
previous measurements and also compatible with the hypothesis of no CPV. They 
bring down the world average for y C p from (1.064±0.209)% to (0.866±0.155)% (7), 
relaxing the previously existing tendency for yep > y [3j . The new world average 
for A r is (-0.022 ± 0.161)% 0. 



Table 1. Recent experimental results for the quantities yep and Ar- The h~h+ yield is given for 
K~ K + and inside the parentheses for tt~it + . First uncertainty is statistical, second is systemat- 
ical. BaBar measures, instead of Ar, AV = (1 + ycp)Ar- 



Experiment 


K~7V + 


h~h+ 


yep (%) 


A r (AY) (%) 




yield (xlO 6 ) 


yield (xlO 3 ) 






LHCb 


0.225 


30 


5.5 ±6.3 ±4.1 


-5.9± 5.9± 2.1 


BaBar 


7.312 


633 (65) 


0.72 ±0.18 ±0.12 


0.09 ±0.26 ±0.09 


Belle 


2.61 


242 (114) 


1.11 ±0.22 ±0.11 


-0.03 ±0.20 ±0.08 



3 Time integrated CP asymmetries 

3.1 Two-body decays 

One of the most discussed results from the LHC in the past year was the first evi- 
dence of CPV in the charm sector, provided by the difference of the time integrated 
CP asymmetries measured in the decays of the D° meson to K~K + and 7r _ 7r + , 



presented by LHCb 11 and supported by new results by CDF 20 and Belle 21 
To first order of approximation, the raw asymmetry can be written as a sum of 
small terms 

A raw (f) = A CP {f) ± A det (f) + A det (n s ) + A prod (D*+), (5) 

where Acp(f) is the physical CP asymmetry of interest, Ad e t(f) is the asymmetry 
in the detection efficiency for the decay to final state /, Ad e t(^s) is the recon- 
struction efficiency for the slow pion from the D* + —> 7r+_D° tagging chain and 
A pro d{D* + ) is a physical asymmetry for D* + production. This asymmetry van- 
ishes in CP conserving strong collisions, as the pp collisions at the Tevatron, if the 
events are reconstructed symmetrically in the pseudo-rapidity r\. However, it is not 
expected to be null in the pp collisions at the LHC. In e + e~ colliders, there is a 
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forward-backward asymmetry in the production of charm mesons arising from jZ 
interference and higher order QED effects [22| . 

When dealing with two-body decays of a spin particle to self-conjugate final 
states /, the detection efficiency Ad e t{f) is null. Since the remaining terms, apart 
from Acp(f), are independent of the final state, they cancel out in the difference of 
asymmetries AA raw — A raw (K + K~) — A ralu (7r + 7r~), which thus becomes a robust 
measurement of the difference between the physical CP asymmetries: 

AA CP = A CP (K+K-) - A cp (tt+tt-) = AA raw . (6) 

Experimentally, the raw asymmetries arc obtained from 

N(D^f)-N(D^f) 
raw[J> - N(D°^f) + N(D°^f) { > 

where N(X) corresponds to the number of reconstructed decays X after background 
subtraction. 

LHCb performed this measurement using 0.62 fb _1 of data collected in 2011. 
Maximum likelihood fits to the dm — m(-K s hr h + ) — m(h~hr) spectra provide the 
number of signal events used to calculate the raw asymmetries in different kine- 
matic bins, in order to guarantee that the detection and production asymmetries 
are cancelled. The kinematic bins are defined in intervals of the -D candidate 
transverse momentum pt and r\ and the final result, shown in Table [2j is obtained 
from a weighted average over this kinematic space. This result deviates from zero 
with a significance of 3.5 a. 

CDF also performed a new AAcp measurement, using the full data sample 



collected at Tevatron Run II, which corresponds to 9.7 fb 1 20 . Compared to the 



previously published measurement 23 , their new result relies on a larger data sam- 
ple and an optimized selection procedure. In order to take into account differences 
in the instrumental asymmetries due to phase space variations among the two final 
states, a reweighting of the samples is performed, so as to equalize the n s impact 
parameter, p? and r/ distributions. A simultaneous % 2 fit of the resulting D°tt + 
and D°ir~ mass distributions provides the number of signal events used in the raw 
asymmetries. The result, also shown in Table [2] is 2.7a away from zero, with a 
precision level very close to the obtained by LHCb. The individual asymmetries, 
measured with a sample of 6.0 fb -1 , are A C p(K-K+) = -0.24 ±0.22 ±0.09% and 



23 



Acp(7r-7r+) = (+0.22 ± 0.24 ± 0.11)% 

In July 2012, Belle presented a preliminary update of the individual asymme- 
tries Acp(K~ K + ) and Acp (7r~7r + ), reporting also a AAqp measurement, using a 
data sample corresponding to 976 fb -1 . Instrumental asymmetries are obtained 
using control samples of CP conserving decays and the production asymmetry 
is eliminated by making the measurement in bins of cosO* , where 9* is the po- 
lar angle of the D* + production in the center of mass system. The results are 
A CP {K-K+) = -0.3 ± 0.21 ± 0.09% and Acp(tt-tt+) = (+0.55 ± 0.36 ± 0.09)%. 
The AAcp value, given in Table [2j deviates from zero with a 2.1er significance. 

The asymmetry difference can be written, to first order of approximation, in 
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Table 2. Recent experimental results for AAqp. The last column shows the difference between 
the average proper time of the K~K + and tt~tt + samples used to extract AAcp- 



Experiment 


K~K+ 
yield (xlO 6 ) 


TV 7T + 

yield (xlO 6 ) 


SA CP 


A(*> 


LHCb 


1.44 


0.38 


(0.82 ±0.21 ±0.11)% 


~ 0.098r D o 


CDF 


0.55 


1.2 


(-0.62 ±0.21 ±0.10)% 


~ 0.26r D o 


Belle 


0.282 


0.123 


(-0.87 ±0.41 ±0.06)% 







Figure 2. Combined measurements of AA^p and Ay, where the bands represent ±1ct intervals. 
The point of no CP violation (0, 0) is shown as a filled circle, and two-dimensional 68% CL. 95% 
CL, and 99.7% CL regions are plotted as ellipses with the best fit value as a cross indicating the 
one-dimensional uncertainties in their center. Reproduced from [3|12|. 



terms of the direct and indirect components as 123] 



AA CP « Aa d c " P + ^-a md , (8) 

T D 

where (t) is the average reconstructed decay time, a-cpif) = (IA/1 2 — 
\A f \ 2 )/(\A f \ 2 + \A f \ 2 ) and a™* = -\ [(|*| - |||) j/coj^- (||| ± |f |)x S in^] and 
the phase <f> is assumed to be universal. From this equation and the A(i) /rrjo values 
given in Tab.[2]for the different experiments, it becomes clear that AAqp is mostly 
a measurement of direct CPV. 

Including higher order corrections, the first term in the right side of Eq. [8 1 24 



changes to Aafj r P [l + ycp{{t{K~ K+)) ± (t(R- K+)))/(2t d o)]. Assuming that the 
current precision of A-p is such that it is not sensitive to Aafjp, HFAG fits the 
contributions of direct and indirect CPV to the most recent results of yep, and 
AA CP by Babar, Belle, CDF and LHCb |l2:, as shown in Fig. [2] This fit yields 
Aafj P = (-6.78 ± 1.47) x 1CT 3 and = (0.27 ± 1.63) x 10" 3 , again favouring a 
direct CPV effect over an indirect, and is consistent with CP conservation with a 
confidence level of only 0.002%. 

BaBar and Belle searched for CPV in the decays of the D + and Df mesons 
to tw o-body final states containing a Kg using data samples corresponding to 469 
fb _1 25|27 and 673 fb _1 26 , respectively. The results are presented in Table [3] 



The first two decays are SCS decays, while the latter result from the coherent 
sum of Cabibbo favoured and doubly Cabibbo suppressed transitions and are not 
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expected to exhibit CPV in the charm dynamics. For all those decays, however, the 
SM predicts a small CP asymmetry due to the K° — K° mixing. This contribution is 
measured to be (+0.332±0.006)% for the D+ -> K° s ir + mode and (-0.332+0.006)% 
for the other modes 29 . The Belle result for the D + — > K®ir + asymmetry was 
performed with a larger sample, corresponding to the full data set of 977 fb _1 , 
and deviates from zero with a significance of 3.2o\ In this work, they estimate 
a correction to the K° — K° mixing induced asymmetry following 30 and, after 
subtracting this contribution, they estimate a CP asymmetry due to the charm 
dynamics alone of (— 0.024±0.094±0.067)%, which is the most precise measurement 
of Acp in charm decays to date. The Acp values obtained for the other channels 
are all consistent with zero and with the SM expectation. In particular the Babar 
measurements presented in May 2011 are the most precise 
these asymmetries. 



27 determination of 



Table 3. Most recent experimental results for CP asymmetries in the decays D+ — > KgK + , 
D+ -> K°n+, D+ -¥ K°n+ and D+ -> K°K+. 



Decay BaBar A C p{%) Belle A C p(%) 


D+ -> K" S K + +0.13 ±0.36 ±0.25 
D+ -> K%-k + +0.6 ±2.0 ±0.3 2 
Df -> K° S K + -0.05 + 0.23 + 0.24 
D + -> K° s -k + -0.44 + 0.13 + 0.10 


27 

7 

27 
25 


-0.16 + 0.58 + 0.25(26 
+5.45 + 2.50 ±0.33|| 
+0.12 + 0.36 + 0.22^ 
-0.363 + 0.094 + 0.067: 
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3.2 Multi-body decays 

In the decays of D mesons to final states containing more than 2 hadrons, the rich 
resonant structure of the intermediate state provide the strong phases necessary 
to the observation of CPV, under the hypothesis of existing weak CPV phases. 
The interference between the resonances can lead to CPV asymmetries which vary 
across the phase space. These local asymmetries can be large compared to the 
overall asymmetry or can even happen in the absence of phase space integrated 
asymmetry. 

LHCb used a model independent approach to look for a CPV signal in the 
SCS decay D + —> K~K + ir + using a high purity sample (~90%) containing about 
370xl0 3 decay s selec ted from about 35 pb -1 of data collected in 2010 |3l] . Within 
this technique 32|33 , the Dalitz plot is divided in bins and for each bin the statis- 



tical significance of the asymmetry 

Nl - aN l 

S 1 cp = , + (9) 
JN l + + a 2 Nl 

is calculated, where N± is the number of D decays in each bin i and a is the 
ratio between the total number of D + and D~ events and is used as a correction 
due to a global production asymmetry. The distribution of S l CP is normal under 
the hypothesis of CP conservation. A y 2 test using x 2 — YliLx &cp provides 
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a numerical evaluation for the degree of confidence for the assumption that the 
differences between the D + and D~ Dalitz plots are driven only by statistical 
fluctuations. Four different binning schemes are considered, two of them uniform 
and two adaptative which account for the resonant structure of the decay. Higher 
statistics control samples consisting of the CF decays Df — > K~K + n + and D + — > 
K~Tt + Tr + are used to demonstrate that possible instrumental asymmetries varying 
over the Dalitz plot are under the sensitivity of the method. The method is also 
applied to the events in the side bands of the D + mass distribution to prove that 
no effect can be introduced by the background. The p-values obtained with the \ 2 
tests for the S % CP distributions using the four different binning choices are 12.7%, 
10.6%, 82.1% and 60.5%. The widths and means obtained from the fits to the 
distributions are all consistent with 1 and 0. No evidence for CPV is found in any 
of the binning schemes. 

BaBar has also performed a similar analysis of D + — > K~K + ir + events 
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using a sample of 476 fb 1 of data with ~ 228 x 10 3 events (92% pure). They 
divide the Dalitz plane into 100 equally populated bins and find a 72% probabil- 
ity for the consistency with the assumption of no CPV. The measured integrated 
asymmetry is A C p — (0.35 ±0.30 ±0.15)%. A full amplitude analysis provides CP- 
violating differences for the phases and magnitudes of the intermediate resonances, 
all consistent with zero. 

CDF searched for CPV in the resonant substructure of D* + tagged D° 
Kgir~ir + decays using a sample corresponding to 6 fb -1 collected at Run 11 35 



The SM predictions for this channel are of the order of 10~ 6 , dominated by the 
Kg mixing contribution. The selected sample is about 90% pure and contains 
approximately 350 x 10 3 events. A simultaneous maximum likelihood fit to the D° 
and D° Dalitz plots, taking into account the effect of efficiency over the Dalitz space, 
provide resonances parameters and CP-violating fractions, amplitudes and phases. 
Differences between the px distributions of the slow pions of opposite charges are 
used to reweight the D° distribution. All the asymmetries are measured to be 
consistent with zero, as is the overall integrated asymmetry Acp = (0.05 ± 0.57 ± 
0.54)%. A model independent analysis similar to the one described above for the 
D + — > K~ K + n + decay also provides null CPV signal. 

For the decay D° — > ir~ n + ir~ ir + , the phase space becomes 5-dimensional. LHCb 
performed the first model independent CPV analysis of this decay, using about 
180 x 10 3 D* + tagged candidates, selected from a sample of 1 fb -1 collected in 
The decay D° — > K~ir + Tr~ir + is used as control mode and events in the 
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2011 _ 

D° sample are rejected in a way such to ensure that the rj and pr distributions of 
the D° decays match those of the D° decays. The distribution of the asymmetry 
significances across the phase space is shown to be compatible with the hypothesis 
of CP conservation. 



4 Summary and Conclusions 

Thanks for the pioneering work of the charm factories and the excellent work done 
by Babar, Belle and CDF, mixing in the D° — D° system is well established, though 
no single measurement with a significance above 5 standard deviations was available 
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until the beginning of the preparation of these proceedings. In November 2012, 
LHCb has submitted a paper reporting the observation of oscillations in the time 
dependent rate of D° — > K + ir~ relative to D° — > K~ir + , which excludes the no 



mixing hypothesis with a significance of 9.1 tr 37 



The recent results obtained for AAcp raised the interest about CPV in charm 
decays. A lot of effort was put, from the theory side, to understand on whether 
this can be, if confirmed, due to a new physics effect or to an enhancement of SM 
penguin contributions [38| . Until now, no clear picture emerged yet. More precise 
measurements of the difference of asymmetries and of the individual asymmetries 
Acp{K~ K + ) and Acp{ir~Tr + ) are needed. Advances in the theoretical under- 
standing of charm decays can also help to disentangle this puzzle, as well as precise 
measurements in additional neutral and charged decay modes. CDF, BaBar and 
Belle ended operation and are finishing to analyse their samples. LHCb has shown 
the ability to reconstruct charm decays with purities comparable to the experiments 
at e + e~ machines and has already collected the largest charm sample of the world, 
corresponding to a luminosity of 1 fb _1 at y/s = 7 TeV and 2 fb _1 at = 8 TeV. 
Precisions comparable or below the SM predictions will only be achievable at the 
LHCb upgrade 39 and the super B factories 40|41 . 
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